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ABSTRACT 

We give further considerations on the problem of the evolution of a coronal, force-free magnetic field 
which threads a differentially rotating, conducting Keplerian disk, extending the work of Li et al. (2001). 
This situation is described by the force- free Grad-Shafranov (GS) equation for the flux function \E'(r, z) 
fSJ ' which labels the poloidal field lines (in cylindrical coordinates). The GS equation involves a function 

, H{'^) describing the distribution of poloidal current which is determined by the differential rotation 

' or twist of the disk which increases linearly with time. We numerically solve the GS equation in a 

sequence of volumes of increasing size corresponding to the expansion of the outer perfectly conducting 
boundaries at (i?m, Zm)- The outer boundaries model the influence of an external non- magnetized 
plasma. The sequence of GS solutions provides a model for the dynamical evolution of the magnetic 
field in response to (1) the increasing twist of the disk and (2) the pressure of external plasma. We 
ly-^ . find solutions with magnetically collimated Foynting jets where there is a continuous outflow of energy, 

' angular momentum, and toroidal magnetic flux from the disk into the external space. This behavior 

contradicts the commonly accepted "theorem" of Solar plasma physics that the motion of the footpoints 
of a magnetic loop structure leads to a stationary magnetic field configuration with zero power, angular 
momentum, and flux outflows. 

In addition we discuss magnetohydrodynamic (MHD) simulations which show quasi-stationary col- 
' limated Foynting jets similar to our Grad-Shafranov solutions. In contrast with the Grad-Shafranov 

solutions, the simulations show a steady uncoUimated hydromagnetic (non force-free) outflow from the 
outer part of the disk. The Foynting jets are of interest for the understanding of the jets from active 
' galactic nuclei, microquasars, and possibly gamma ray burst sources. 

o ■ 

' 1. INTRODUCTION equation are calculated inside given conducting bound- 

r\ ' tj- 1,1 IT 4- J -4- 1 J- J- J • + u J aries. The solutions are self-consistent in the respect that 

tii Highly-collimated, oppositely directed lets are observed ,, ; . , r i r. i i i- • ;i ; i ; ;i tKt .■ i 

I ' ■ °. , . J / r 1 -ji (1 the twist of each held line is that due to the diflerential 

Q , m active galaxies and quasars (see tor example Lsridle & otat'o of a Ke le 'a d'sk 

Eilek 1984), and in old compact stars in binaries (Mirabel t nr. • / i i • i- 

S & Rodriguez 1994; Eikenberry ct al. 1998). Further, weU J'^Ji we summarize the theory, and m §3 we discuss 

^ ' 11- .1 • • r • 4. • J- 11 1, the ditterent types ot solutions found. In b4 we develop 

rrt . collimated emission line lets are seen m young stellar ob- ,,. , , r .■ ■ . xpr i- ,i 

^ ■ jects (Mundt 1985; Biihrke, Mundt, & Ray 1988). Recent ^^^^^^^^ ^"f^^ Poy^t/ng jets. In §5, we discuss the 

K>. ,V , -f J J t-- consequences ot expanding the boundaries. Sections 6-9 

work tavors models where twisting 01 an ordered magnetic ,. , f , tpi« 

k j2 1 J ii J- J.- J- 1 J. J. 11 discuss the relevant conservation laws. In qlU we compare 

. held threading an accretion disk acts to magnetically ac- ,, t-, ,. ,n -.i , •<■ n i i i • i 

1^ ■ celerate the jets (Meier, Koide, & Uchida 2001). There he Poyntmg outflows with centrifugally launched winds. 

rH • . . /iN.u Lj u Inhllwe discuss the collapse ot the inner part ot the disk 

'■^ , are two regimes: (1 the huaromaqnetic reqime, where en- , , ,n n -- 

■ " " J 1 i ■ J u u 1 due to angular momentum outflow to Foynting lets. In 

ergy and angular momentum are carried by both the elec- ... „ r t-. • 

, c ij J ii. 1 • i- £1 f 1 ■ 1 ■ bl2 we discuss conditions tor occurrence ot Foynting lets 

tromagnetic field and the kinetic flux ot matter, which is . , ,. ,, . n <• ,i i • i • . i \ no 

r ,11 u • -L J /r>\ including the influence oi the kmk mstabiiity. In we 

relevant to the lets irom young stellar objects; and (2) . °, tutttt^ ■ i ■• i • i • ti , 

,x. n a J If give results oi MHD simulations which give Foynting lets 

the Foyntmq flux reqime. where energy and angular from . .i ^ i m r i n , . ■ 

,1 J. f • J J • ii u 1 i- similar to our Grad-Shairanov solutions. Section 14 gives 

the disk are carried predominantly by the electromagnetic i . r , i . i 

cij i,-i,- 1 i li- J • conclusions from this work, 

neid, which is relevant to extra-gaiactic and microquasar 

jets, and possibly to gamma ray burst sources. 

Stationary Foynting flux dominated jets have been 2. THEORY OF POYNTING OUTFLOWS 

found in axisymmetric MHD simulations of the opening Here, we consider further the theory of Foynting out- 

of magnetic loops threading a Keplerian disk (Romanova flows (Li et al. 2001). We assume that magnetic fleld 

et al. (1998); Ustyugova et al. 2000). Theoretical studies loops threads a differentially rotating highly conducting 

have developed models for Foynting outflows from accre- Keplerian accretion disk at some initial time t = 0. Above 

tion disks (Lovelace, Wang, & Sulkanen 1987; Colgate & Li the disk we assume a "coronal" or force-free magnetic fleld 

1998). The present work represents a continuation of the in the non-relativistic hmit. This situation is described by 

study by Li et al. (2001) where self-consistent, axisym- the force-free Grad-Shafranov (GS) equation for the flux 

metric, non-relativistic solutions of the Grad-Shafranov function ^{r,z) which labels the poloidal field fines (in 
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Fig. 1. — "Initial" dipole-like vacuum magnetic field. In this and subsequent plots, r and z are measured in units of the 
radius ro of the O— point in the disk plane. The solid lines are the magnetic field lines for the case where the flux function 

on the disk surface is 5'(r. z = 0) ~ c^r^Kjia^ + r^)^/^ with a = ro/\/2. In this and subsequent plots 'I' is measured in 

units of *o = rlK/B^/"^. Note that B^{0) = 6V3*o/ro « 10.4*o/ro. The dashed lines are the field lines for the case where 
the outer boundaries {Rm = 10, Zm = 10) are perfectly conducting; for this case ^'(r, z = 0) ^ ^{r, 0)[1 — (r/ro — 1)^/81] 
so that the O— point is still at ro and is unchanged. Because of axisymmetry and reflection symmetry about the z = 
plane, the field need be shown in only one quadrant. 



cylindrical coordinates). The GS equation involves a func- 
tion i?(^') describing the distribution of poloidal current 
which is determined by the differential rotation or twist 
of the disk which increases linearly with time. We nu- 
merically solve the GS equation in a sequence of volumes 
of increasing size corresponding to the expansion of the 
outer perfectly conducting boundaries at (i?,„, Zm). The 
outer boundaries model the influence of an external non- 
magnetized plasma. The sequence of GS solutions provides 
a model for the dynamical evolution of the magnetic field 
in response to (1) the increasing twist of the disk and (2) 
the pressure of external plasma. 

Cylindrical (r, </>, z) coordinates are used and axisym- 
metry is assumed Thus the magnetic field has the form 
B = Bp + B^cj) , with Bp = Brr + B^z . We can 
write Br = -{l/r)(d'i' /dz), B^, = (l/r)(9*/ar), where 
^{r,z) = rAcf,{r,z). In the force- free limit, the magnetic 
energy density B^/(87r) is much larger than the kinetic or 
thermal energy densities; that is, the flow speeds are sub- 
Alfvcnic, v^ <C = B^/47rp, where va is the Alfven ve- 
locity. In this limit, J x B; therefore, J = AB (Gold 
& Hoylc 1960). Because V- J = 0, (B ■ V)A = and 
consequently A = A(^'), as well-known. Thus, Ampere's 
equation becomes V x B = 47rA(\l/)B/c. The r and z 
components of Ampere's equation imply 

rB^ = H{^), and dH{^)/d^ = A'k\{^)/c, 

where H{^) is another function of Thus, H{^) = 

const are lines of constant poloidal current density; Jp = 
(c/47r)(dff/d*)Bp so that (Jp • V)H = 0. The toroidal 
component of Ampere's equation gives 



(1) 



with A* = d'^/dr^ - {l/r){d/dr) + d'^/dz^, which is the 
Grad-Shafranov equation for 4* (see e.g. Lovelace et al. 
1987; Li et al. 2001). 

Ampere's law gives ^dl • B = (47r/c) J dS ■ 3, so 
that rB^{r,z) = -ff(^') is (2/c) times the current flow- 
ing through a circular area of radius r (with normal 
z) labeled by '^{r,z)= const. Equivalently, — Jf[\l/(r, 0)] 
is (2/c)x the current flowing into the area of the disk 
< r. For all cases studied here, —H{^) has a maxi- 
mum so that the total current flowing into the disk for 
r < r„i is Itot = {'2/c){-H)max, where is such that 
—H['^{rm,0)] = {—H)max SO that rm is less than the ra- 
dius of the O— point, tq. The same total current Itot flows 
out of the region of the disk r = to tq . 

We consider an initial value problem where the disk at 
t = is threaded by a dipole-like poloidal magnetic field. 
The form of (4") is then determined by the differential ro- 
tation of the disk: The azimuthal twist of a given field line 



going from an inner footpoint at r\ to an outer footpoint at 

r2 is fixed by the differential rotation of the disk. The field 
line slippage speed through the disk due to the disk's finite 
magnetic diffusivity is estimated to be negligible compared 
with the Keplerian velocity. For a given field line we have 
rd0/_B0 = dsp/Bp, where dsp = V dr^ -|- dz'^ is the poloidal 
arc length along the field line, and Bp = -^/B^ + B^. The 
total twist of a field line loop is 



i: 



dSr, 



—B^ 



dsp 



(2) 



with the sign included to give Acj) > 0. For a Keplerian 
disk around an object of mass M the angular rotation is 
ljk = \J GMj r^ so that the field line twist after a time t 
is 



A(j!)(*) =(Jot 



3/2 



3/2' 



(wo t) J^(*/*o) 



where ro is the radius of the O— point, wq = \fGMjr\ 



0' 



(3) 
and 

^ is a dimensionless function (the quantity in the square 
brackets). The O— point is the point in the midplane of 

the disk encircled by the poloidal magnetic field lines; at 
this point Bp = 0. At sufficiently small r\ one reaches the 
inner radius of the disk rj (<i; ro) where we assume lok 
saturates at the value loki = \/ GM /r^. For the dipole- 
like field of Figure I, k. 39/*(*o/*)^/'' for ^/^o < 1, 
while T w 3.64(1 - ^/^of/^ for 1 - */*o < 1- 

For an accretion disk aroimd a massive black hole 
M = Mf,lQP Mq in the nucleus of galaxy, the twist is 

T = _(</3.17d)(ro/10i5cm)3/2//M8- The inner radius of 
the disk is r, ~ MgQ x 10^^ cm for a Schwarzschild black 
hole. 

3. NUMERICAL SOLUTIONS OF GRAD-SHAFRANOV 
EQUATION 

We solve equation (1) inside a cylindrical "box," r = to 
Rm and 2 = to Zm with H{^) determined by iteratively 
solving equations (1) - (3) (see Li et al. 2001; Finn & Chen 
1990). The outer boundaries (r = Rm and z = Zm) are 
treated as conducting surfaces representing the interface 
between the coronal B field and external, non-magnetized 
plasma. The external plasma will be pushed outward by 
the coronal field as it is twisted. Thus we consider the B 
field behavior as Rm oo and Zm — * oo. 

For a dipole-like field threading the disk (Figure 1), it 
is natural to measure lengths in terms of the radius of the 
O— point in the disk, ro. The natural value of the flux func- 
tion ^! is its maximum value at the O— point, . In turn, 
a natural unit for the magnetic field strength is 4*0 /^q. 
Notice that for the case of Figure 1, ^o/^o ~ Bz(0)/10.4. 
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Fig. 2. — Summary of dependences of numerical solutions of equation (1) on the twist T = uiot. For these solutions 
-ff satisfies equations (2) and (3). Here, Wm is the magnetic energy (with WmO ~ 0.57 the energy for zero twist), $t is 
the toroidal flux, and {—H)max is the total poloidal current (times c/2). The solutions are inside a conducting box with 
Rm = 10 = for an initial dipole-like field with O— point at Tq = 1 as shown by the dashed lines in Figure 1. The value 
of T where the curves jump separates the two types of solutions found. The error bar indicates the estimated uncertainty 
in the values at T 2 which results from the finite grid used in calculating A0(^'). 



Further, a natural unit for the total current (times c/2), 
c/tot/2 = {—H)jnax is ^o/rQ. Also, we measure the total 
toroidal flux and the magnetic energy, 

$t = y" drdz B^, Wm = \j rdrdz , (4) 

in units of and ^'o/'^O) respectively. For boundary con- 
ditions, on the disk surface, \l/(r, 0) is fixed and equal to 
its "initial" value on the disk owing to the disk's perfect 
conductivity. On the z— axis wc can take ^(0,z) = 0, 
because B^ = —{l/r){d'^/dz) = 0. On the conducting 
outer boundaries, ^{r,Zm) = and ^{Rm,z) = 0. Note 
however that in §14 we discuss simulation results for the 
case of open outer boundaries. A uniform (r, z) grid of 
200 X 200 was used. 

Figure 2 shows the behavior of a set of field solutions of 
equation (1) as the twist T = uot is increased. The nature 
of the solutions changes dramatically as the twist increases 
above a critical value Tc w 1.14 rad. For T < Tc the na- 
ture of the field solutions is shown in Figures 3 and 4. The 
twisting of the field by the differential rotation of the disk 
"pumps" magnetic flux and energy into the disk corona 
and the field tends to "inflate.'' This behavior of coronal 
magnetic field loops of the Sun as a result of footpoint 
twisting is well known from the works of Aly (1984, 1991) 
and Sturrock (1991). The self-similar inflation of a force- 
free fleld threading a non-Keplerian disk without outer 
boundaries was studied Lynden-Bell and Boily (1994) and 
their solution for small twists is analogous to our low-twist 
solutions. The expansion of a force- free field into flnite 
pressure external plasma has been studied by Lynden-Bell 
(1996) and Li et al. (2001) and the poloidal field is found 
to fairly uniformly fill the coronal space. 

For T > Tc, a new field configuration appears with a 
different topology. This is shown in Figure 5 where it is 
seen that a "plasmoid" consisting of toroidal flux has de- 
tached from the disk and is separated by the dashed field 
line which has an X— point above the O— point on the 
disk. Figure 6 shows a 3D view of two representative field 
lines for the same case. These high- twist equilibria consist 
of a region near the axis which is magnetically collimated 
by the toroidal -B^ field and a region far from the axis, 
on the outer radial boundary, which is anti- collimated in 
the sense that it is pushed against the outer boundary. 
The field lines returning to the disk at r > rg are anti- 
collimated by the pressure of the toroidal magnetic field. 
The poloidal field fills only a small part of the coronal 
space. In a purely analytical analysis, Heyvaerts (2001) 
has independently found MHD equilibria involving the si- 
multaneous formation of a collimated axial jet and an un- 
coUimated outflow. Figure 6 shows a three dimensional 
view of sample fleld lines. 

As a test of our numerical solutions note that conserva- 



tion of axial momentum can be written as 

V • (T • z) = , T i = B^z/Stt - B^'Bp/A-n . (5) 
Integration of equation (5) over the "box" [Rm, Zm) gives 

rdr {Bl + Bl-Bl),=o = rdr {BI),=z^ , (6) 
Jo Jo 

where the integral on the left hand side represents the flux 
of momentum from the disk into the box and the other 
integral the flux of momentum out of the top of the box. 
This equation is accurately satisfled by our numerical so- 
lutions; the typical errors are < 0.1%. 

4. ANALYTIC SOLUTION FOR POYNTING JET 

Most of the twist Acf) of a field line of a Poynting 
jet occurs along the jet from z = to Z^- Because 
-r'^d(j)/H{^) = dz/B^, wc have 

A0(^) ^ {^^t)T{^/^o) _ Z^ 

-H{^) -F(«') ~ r'^B, ' ^ ' 

where r'^Bzir,z) is evaluated on the straight part of the 
jet at r = r('I'). In the core of the jet <C V&o, we 
have « 3^/*(\E'o/^')^^'^, and in this region we can take 
* = C'I'o(r/ro)«, and H = -/C(*o/^'o)(*/^'o)^ where 
C, q, JC, and s are dimensionless constants. Equation 
(1) for the straight part of the jet implies q = 1/(1 — s) 
and C2(i-«) = s(l - s)2/CV(l - 2s). Thus we find 
s = 1/4 so that q = 4/3, C = [9/32]2/3/c4/3^ and 

IC = 3yH{rooJot/Zm). 
In order to have a Poynting jet, we find that /C must 

be larger than 0.5. This corresponds directly to the 
condition T > Tc for the occurrence of the high-twist so- 
lutions. The condition arises from the fact that there is 
a competition between the build up of toroidal flux inside 
the box due to twisting by the disk which acts to increase 
B(j, and the expansion of the boundaries which acts to de- 
crease i?0 (sec §12). If the boimdaries expand too rapidly 
B(j, does not increase sufficiently to give a self-collimated 
Poynting jet. For the case of uniform expansion of the 
top boundary, Zm = V^t, this condition is the same as 
V, < 9.2(roWo)- For the case of Figure (5), /C w 0.844. 
The fleld components in the straight part of the jet are 




These dependences agree approximately with those found 
in numerical simulations of Poynting jets (Ustyugova et al. 
2000). On the disk, * w 33/2*o(r/ro)2 for r < ro/33/4. 
Using this and the formula for \l/(r) gives the relation be- 
tween the radius of a field line in the disk, denoted r^. 
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Fig. 3. — Poloidal field solutions for case of conducting boundaries for twists T = 0.67 rad. (dashed lines) and T = 1.1 
rad. (solid lines). The initial poloidal magnetic field is shown by the dashed lines in Figure 1. 



Fig. 4. — The top panel shows the field line twist function ^{^), with the circles indicating the derived numerical values 
and the smooth curve the theoretical dependence for a Keplerian disk given by equation (3). The case shown corresponds 
to the solid field lines in Figure 2 where the twist is T = 1.1 rad. Here, ^ is measured in units of ^o, the value at the 
O— point. For this case, the lower panel shows the poloidal current flow — _ff(^') (in units of '^o/ro), with the points 
indicating the numerical values and the smooth curve an analytic fit. The initial poloidal magnetic field is shown by the 
dashed lines in Figure 1. 



and its radius in the jet, r/ro = 6.5{rd/rQ)^^'^IC~^ . Thus 
the power law for 4* is applicable for ri < r < r2, where 
ri = 6.5ro(ri/ro)'^/^//C and r2 = 1.9ro//C, with n the in- 
ner radius of the disk. The outer edge of the Poynting jet 
has a transition layer where the axial field changes from 
Bz{r2) to zero while (minus) the toroidal field increases 
from —B^{r2) to {—H)max/f2- Using equation (6), which 
is only approximate at r2, gives {—H)max ~ l-2/C5'o/ro, 
This dependence agrees approximately with our Grad- 
Shafranov solutions. 

5. EXPANSION OF BOUNDARIES 

The magnetic forces on the outer wall increase by a large 
factor in going from the low- twist to high- twist solutions. 
The radial force on the cylindrical wall and the axial force 
on the top wall are 



1 



dz Bi 



rdr Bl 



(9) 

For the low- twist solution of Figure 3 (T = 1.1), (i^V, J^) w 
(0.0061,0.013), whereas for the high-twist solution of Fig- 
ure 5 (T = 1.79), {Fr,F^) w (0.26,0.45) in units of 

(^o/ro)^ 

Figure 7 shows the radial dependence of the magnetic 
pressure on the top boundary for the cases of low-twist 
and high- twist solutions. 

Different behavior is exhibited by the low-twist and 
high-twist solutions as the conducting boundaries are 
moved outward. For the low-twist solutions, the poloidal 
field tends to expand outward to fairly uniformly fill the 
available space. As Rm — > oo and Zm — > oo, we find that 
these solutions are similar to those obtained by Lynden- 
Bell and Boily (1994) where there are no outer boundaries. 

In contrast, for the high-twist solutions, the poloidal 
field near the axis maintains its coUimation due to B^, 
whereas the "return" poloidal field far from the axis is 
pushed against the outer cylindrical wall due to the B^ 
field. As Rm ^ oo and Zm oo, the coUimatcd field 
near the axis will extent outward along the z— axis in the 
absence of instabilities. It is clear from Figure 7 that the 
magnetic pressure on the top boundary peaks near the 
axis. Thus, this region of the boundary should expand 
most rapidly in the physical case where the boundary is 
an interface with external plasma. 

Estimation of the axial expansion of a Poynting jet 
can readily be made assuming a region of radius ^ gr2 
(g ~ 2 — 3) of the jet expands with velocity Vz into a 
constant density external medium. For likely conditions 
Vz is much larger than the sound speed in the external 



medium so that the ram pressure due to the jet motion 
is pextyzi assuming (V^/c)^ <C 1 as required by our non- 
relativistic treatment. Balancing this pressure with the 
magnetic pressure of the jet gives 



PextV,^ 



i-H) 



2 

max 



2„2 



Q.U[Bz{Q)/g]\rouJo/Vz, 



or 



Vz^{r,u:,f/^Bz{0)lgY'VpTt 



(10) 



The condition for a Poynting jet /C > 0.5 corresponds to 
Vz < 9.2(roWo)- For B^(0) = lOOG, ro = lO^^cm, and 
M = IO^Mq, the external density Pext must be larger 
than ~ 2 X 10~^^g/cm^ in order to have Vz < c. For larger 
magnetic fields, a much larger external density is needed 
to give Vz < c. This points up the need for a relativistic 
treatment of the Poynting jet expansion. 

6. ANGULAR MOMENTUM CONSERVATION 

Conservation of angular momentum about the 2;— axis 
can be written as 

V • J = , J = -rBpB^/4TT = -BpH/4TT , (11) 

where J is the angular momentum flux-density vector. In- 
tegration of equation (11) over the "box" {Rm, Zm) gives 



'■R„ 



rdr H{Bz)z 



=0 • 



(12) 



The subscript z = here and subsequently indicates 
that the quantity is evaluated on the top surface of the 
disk. For a dipole-like field threading the disk where 
^-(^^,0) = = *(0,0), equation (12) gives 



-- d^H{^)-- d^H{^) . 



(13) 



The flrst integral represents the outflow of angular mo- 
mentum from the inner part of the disk (interior to the 
O— point), and this equals the angular momentum inflow 
into the outer part of the disk given by the second integral. 

The outflow of angular momentum from the inner part 
of the disk causes enhanced accretion of this part of the 
disk, whereas the inflow of angular momentum to the outer 
disk reduces the accretion rate. Because the Poynting out- 
flows carry negligible matter, the continuity equation for 
the disk is 

dE/dt + \7iT,Vrr) =0 , (14) 



5 



Fig. 5. — Poloidal field lines for Poynting jet case for twist T = 1.84 rad. and {—H)max = 1-13. The initial poloidal 
magnetic field is shown by the dashed lines in Figure 1. The dashed contour is the separatrix with the X-point indicated. 
Note that the radial width of the upgoing field lines along the axis is about one-half the width of downgoing field lines at 
the outer wall as required for equilibrium. 

Fig. 6. — Three dimensional view of two field lines originating from the disk at a; = ±0.32 ro = 0.4) for the Poynting 
jet of Figure 5. Each field line has a twist of 8.22 rad. or about 1.31 rotations about the 2;— axis from its beginning at 
ri and end at r2- The z— axis is tilted towards the viewer by 30°. 

Fig. 7. — Radial dependences of the magnetic pressure Pmag = {r, Z^) / 4:TT on the top boundary for the cases of 
high-twist (T — 1.79) (the top curve) and low-twist (T = 1.1) (bottom curve) solutions. The pressure is in units of 
i^o/rir « (S.(0)/10.4)2. 



where T,{r,t) is the surface mass density of the disk. The 
continuity equation for the disk angular momentum is 

d{Ee)/dt + V-{Evre r + T;'^' f) = r{B^B,),=o/2-K , (15) 

where 7^"'* represents the viscous transport of angular mo- 
mentum in the disk, the term on the right-hand side the 
outflow or inflow of angular momentum (from the two sides 
of the disk) due to the magnetic field, and £ is the specific 
angular momentum of the disk matter. The disk is as- 
sumed almost Keplerian so that i = \/GMr = rvx and 
consequently the two continuity equations give the mass 
accretion rate M — Mb + Mms , where 

MB{r) = -2nrEvr = -2{r^ Ivk){B^B,),=o , (16) 

is the "magnetically driven" mass accretion rate, and 
Myis = 67ri/r d{i'^/rT,)/dr is the viscous accretion rate 
with v the kinematic viscosity (Lovelace, Newman, & Ro- 
manova 1997). We have Mb > (or< 0) for r < ro (or 
r > ro). The accretion speed is u = —Vr = ub + Ums with 
Ub = MB/(27rrS) and u^,s = Mvts / i^Tir^) . ^ 

Figure 8 shows the radial dependence of Mb for a high- 
twist case. That Mb due to the Poynting outflow is a 
function of r emphasizes the fact that disk is not station- 
ary. 

7. MAGNETIC FIELD TRANSPORT IN THE DISK 

Poloidal magnetic field threading the disk tends to be 
advected inward with the accretion flow, but at the same 
time it may diffuse through the disk owing to a finite mag- 
netic diffusivity rjm of the disk. The continuity equation 
for poloidal flux through the disk Bz{r, 0,t) is 

dB,/dt + V-[vrB,r + UrB.r] = , (17) 

where Ur — {rjm/h) tan(0) is the outward diffusive 
drift speed, h the half-thickness of the disk, tan(^) = 
{Br/Bz)z=o, and a smaller, second order diffusion term 
iVrnd'^Bz/dr'^) has been omitted (see Lovelace et al. 1997). 
For cases where the diffusivity is of the order of the vis- 
cosity and where the viscosity is given the Shakura and 
Sunyaev (1973) prescription v = acsh (with a < 1 and 
Cs the midplane sound speed), the diffusive drift speed is 
Ur ~ acs tan{9). For a strong magnetic fleld threading the 
disk, the accretion speed u in the inner part of the disk 
(r < ro) may be large with ub ^ and ub > Ur so 



that the disk flow advects the ^^(r, 0) field inward. On 
the other hand at large radii, the magnitude of u is prob- 
ably much smaller so that the B^ (r, 0) field tends to drift 
outward, Ur > u. 

8. ENERGY GONSERVATION 

We assume the coronal plasma is perfectly conducting 
so that E = —V X B/c. For quasi- stationary and axisym- 
metric conditions, V x E = 0, and thus = 0,and 
Ep = — V$. so that Vp a ±Bp, and the electro- 
static potential $ = $(4'). Thus E = -fiV^/c with 
f7(5') = d^{'^)/d'^. For the situations considered here, 
all field lines pass through the disk. At the disk surface, 
Er(r,0) = -{v^)diskBz{r,0)/c since = 0. Therefore, 
n[^ir,0)] = {v^)d^sk/r. 

For a force-free plasma we have 




Integration of equation (18) over the "box" (r = to 
Rm, -2 = to Zm) gives 

dWm 1 ^ r R R ^ 

~dr " ~2 X \'"'t'B4>Bz)z=a 
1 . /•^'" 2 

— -rRmaxRm / dz H (Rm, z) 

4 Jo 

--;Zm rdrB'^{r,Zm) , (19) 

4 Jo 

where Wm is the magnetic field energy in the box (equation 
4). The second and third integrals in equation (19) rep- 
resent the energy expended in "pushing" the boundaries 
outward (the pdV work). The first integral represents the 
outflow of energy from the inner part of the disk surface 
(inside the O— point) and the inflow into the outer part of 
the disk. For a dipole-like field threading the disk where 
*(i?„,0) = = ^(OjO), the first integral for the power 
output from the disk can be rewritten as 

-i r°rf*A<?i (^-)F(*)=^^ f d^T[-H{n, 

^ Jo ^'^o Jo 

(20) 

where wo = \/ GM/tq, and the tildes indicate dimension- 
less variables. For the low-twist solution (T = 1.1) of 
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Fig. 8. — The top panel shows the radial dependence of the mass accretion rate of the disk Mb due to the Poynting 
outflow from the disk, and (0.1 x) the Poynting power outflow per unit radius dEg/dr for quasi-stationary conditions. 
The lower panel shows the radial dependence of the angle between poloidal field lines and the z— axis at the disk surface, 
6 = t&n~^ {\Br/Bz\). Both panels are for the high-twist case of Figures 5-7. As mentioned, for ^ > 0.95 or 0.764 < r < 1.33 
the field lines are closed. In the interval 1 < r < 1.5, Br/B^ is less than zero. 



Figure 3, the dimensionless integral denoted X » 1.75, 
whereas for the high-twist solution of Figure 5 (T = 1.84) 
I « 3.05. Thus, E « 



/ ro ( M \ 

VlQiScm/ \10»MqJ 



1/2 



(21) 

for the power output from the two sides of the disk. 

The Poynting power outfiow per unit radius from the 
two sides of the disk is 

dEe/dr = -tok(S^S,)^=o = «/2r) Mb , (22) 

for quasi-stationary conditions. Figure 8 shows the radial 

dependence of dEB/dr which indicates that most of the 
power outfiow occurs in the inner part of the disk. 

9. GENERATION OF TOROIDAL FLUX 

Perfect conductivity and Faraday's law imply 

dB 



dt 



= V X (v X B) 



(23) 



where v is the plasma flow velocity. We apply this equa- 
tion to a "box" extending from the disk at z = to a 
non-rotating, perfectly conducting surface at z = Zm, and 
from the axis r = to a cylindrical, non-rotating, per- 
fectly conducting surface at r = Rm- These outer surfaces 
are allowed to move so that Rm and are in general 
time-dependent. The toroidal flux $t (equation 4) obeys 



d^^f"' 
dt Jo 



dr{v^Bz-VzB^) 



dziVrB^-VijyBr) 



(24) 

On the top boundary z = Zm, B^ = and i?^ = 0; on the 
outer cylindrical boundary r = Rm, B^ = and B^ = 0; 
on the 2;— axis B^ = and Bj. = 0; and on the disk 
Vz{r,z = 0) = 0. Consequently, equation (24) simplifies 
to ^ 

I dr v^Bz{r,z = 0) , (25) 
/o 



dt 



where = ^jGMj r is the azimuthal velocity of the 
disk. The integrand of equation (25) is independent of 
time in view of our assumption that ^'(r, 0), and thus 
Bz{r,Qi) = (l/r)9\E'(r, 0)/9r, is time-independent. 

For a dipole-like field threading the disk [where \I/(r, 0) 
increases from zero to a maximum "^f) and then decrease 
to zero at r = Rm], equation (25) can be rewritten as 



d^t _ /"^° 
dt Jo 



(*) 



(26) 



where (\E') is the time derivative of equation (3), which 
is independent of time. For the dipole-like field of Figure 



1, evaluation of this integral gives d^t/dt « — 12.1a;o^o 
-12.1wo[S^(0)/10.4]rg. Thus, d<^t/dt» 



-1.3 X 10 



28 



Gcm2 f BzjO) 
3kG 



V 1015cm/ 



/2 / M 



V lOSMn 



1/2 



(27) 

is the toroidal flux generation rate from one side of the 
disk. 

10. POYNTING VERSUS CENTRIFUGAL OUTFLOWS 

The centrifugal force near the surface of a Keplerian disk 
has a role in launching hydromagnetic outflows (Bland- 
ford & Payne 1982). The "centrifugal launching" requires 
that the field lines (projected into the poloidal plane) be 
tilted away from the z— axis by an angle 9 greater than 
30°. The magnetic force, dominantly d{B'^/8TT)/dz, is 
comparably important for launching hydromagnetic out- 
flows (Lovelace, Berk, & Contopoulos 1991; Ustyugova et 
al. 1999); it also increases as B increases. Figure 8 shows 
the radial variation of 0{r) = t&n^^ (\Br/ Bz\)z=(h, for a 
high-twist case. For this case, 6 < 30° for r < 0.3ro, which 
includes the part of the disk giving the largest power out- 
put per unit radius, dEB/dr, shown in the top panel of 
Figure 8. Thus there is a Poynting outflow under condi- 
tions where no centrifugal outflow occurs. 

However, for the part of the disk where 9 > 30° we pre- 
dict a centrifugal hydromagnetic outflow with —{B^Bz) > 
and inward magnetically driven accretion, ub > (sec 
Ustyugova et al. 1999). Hydromagnetic outflows from the 
region r > ro are found in the MHD simulations discussed 
by Ustyugova et al. (2000) and in §13. 

11. COLLAPSE OF INNER DISK 

Quasi-stationary Poynting jets from the two sides of the 
disk within ro give an energy outflow per unit radius of 
the disk dEB/dr = rvK{—B0Bz)h, where the h subscript 
indicates evaluation at the top surface of the disk. This 
outflow is ~ rodEB/dr t^A:(r'o)(^'o/?"o)^ is estimated in 
equation (22), which agrees approximately with the values 
derived from the simulations (see §14). 

For long time-scales, the Poynting jet is of course time- 
dependent due to the angular momentum it extracts from 
the inner disk (r < ro). This loss of angular momentum 
leads to a "global magnetic instability" and collapse of 
the inner disk (Lovelace et al. 1997). An approximate 
model of this collapse can be made if the inner disk mass 
Mrf is concentrated near the O— point radius ro{t), if the 
field line slippage through the disk is negligible (see §8), 
^'o =const, and if {—rB^)max ^ ^o/roit) (as found here). 
Then, M^dro/dt = -2^l{GMro)-^'^ ■ If ti denotes the 
time at which ro(fj) = r, (the inner radius of the disk), 
then ro{t) = ri[l - {t - ti)/tcoiiY'^ , for t < ti, where 

tcou = \fGM Mdr,^^V(3*o) is the time-scale for the col- 
lapse of the inner disk. (Note that the time-scale for ro to 
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decrease by a factor of 2is~ ti{ro/ri)^^^ » U for ro ^ rj.) 
The power output to the Poynting jets is 



<3 tcoU 



t-U 

tcoll 



-5/3 



(28) 



where AEioi = GMMd/Iri is the total energy of tlic 
outburst. Roughly, tcoii ~ 2 day MKM^/Mq) (6 x 
10^^Gcm^/\l/o)^ for a Schwarzschild black hole, where va- 
lidity of the analysis requires tcoii ^ ti- Such outbursts 
may explain the flares of active galactic nuclei blazar 
sources (Romanova & Lovelace 1997; Romanova 1999; 
Levinson 1980) and the one-time outbursts of gamma ray 
burst sources (Katz 1997). 



12. 



OCCURRENCE OF POYNTING JETS AND KINK 

INSTABILITY 



The rate at which toroidal flux <I>t is created in the re- 
gion above the disk is d^t/dt « ~\2ujq<^q (sec §10). If the 
area R.„iZ„i of this region is fixed owing to a very dense 
external plasma, the average value of — > in it in- 
creases. On the other hand if the area RmZm increases 
rapidly enough, —B^ will decrease. The Poynting jets oc- 
cur under conditions where —B^ increases to a sufficient 
extent to cause pinching of the poloidal magnetic field. Be- 
cause the most rapid expansion of the B field of occurs in 
the 0— direction, a necessary condition for occurrence of a 
Poynting jet is that the rate of expansion of the boundary 
Vz = dZm/dt be bounded by some constant. In fact the 
condition obtained in §4 for occurrence of Poynting jets 
has this form, Vz < 9.2{LOoro)- 

Note that there may be "self-regulation" in the respect 
that the field configuration which occurs is at the "bound- 
ary" between low and high- twist solutions. In view of Fig- 
ure 7, the low-twist field gives a gradual expansion of the 
boundaries which allows build up of toroidal flux whereas 
the high-twist field gives a more rapid expansion which 
tends to give a slower increase of 

The region of coUimated field (see Figure 5) - the Poynt- 
ing jet - has <C v\ and is kink unstable according 
to the standard non-relativistic analysis (e.g., Bateman 
1980, ch. 6). The instability will lead to a helical distor- 
tion of the jet with the non-linear amplitude of shift of 
the helix Ar ~ VAt and with the helix having the same 
twist about the 2;— ax;is as the axisymmetric B field. Note 
however that for the astrophysical conditions of interest 
VA = |B|/-\/47r/9 is likely to be larger than the speed of 
light. A relativistic perturbation analysis is then required 
including the displacement current. The physical Alfven 
speed is Va = c/(l -|- c^/f^)^/^ < c. Thus the speed of 
lateral displacement of the helix is less than c. The evolu- 
tion of the Poynting jet evidently depends on both Va and 
the velocity of propagation of the "head" of the jet Vz (§5) 
which may be relativistic. Relativistic propagation of the 
jet's head may act to limit the amplitude of helical kink 
distortion of the jet. On the other hand a sub-rclativistic 
propagation of the head may allow the helix amplitude to 
grow but this amplitude can be limited by flux conserva- 
tion as discussed by Kadomtsev (1963). 

13. MHD SIMULATIONS OF POYNTING JETS 

For the MHD simulations described here, the initial 
magnetic fleld has a dipole-like form as shown in Figure 1. 



The computational region r = to i?„ia.r. z = to Z^ax is 
taken to have Rmax = Z^ax ~ IOt'o- Initially, the corona 
of the disk is in isothermal equilibrium without rotation. 
At t = the disk starts to rotate with Keplerian velocity 
V4,{r,0) = rVlK, where VIk = \/GM /(r^ + r^)-''/'*, where 
the smoothing length = 0.2ro is interpreted as the inner 
radius of the disk. The smoothed gravitational potential 
is -GM/(r2 +r2)V2. 

On the disk surface, the boundary conditions are as fol- 
lows (Ustyugova et al. 2000). Two of the boundary condi- 
tions come from the fact that the tangential electric field 
(E')t in the frame rotating with the disk (at the Keplerian 
velocity) is zero; Bz at the disk surface is time-independent 
whereas Br and at the surface vary with time. Two 
further boundary conditions fix the; entropy of the plasma 
coming out of the disk to be Sd(r) and the density of the 
outflowing plasma to be Prf(r). If Vz at the disk surface, 
calculated by solving the MHD equations in the compu- 
tational region, increases to the point where it is larger 
than the slow magnetosonic speed in the z— direction at 
the disk's surface Cgmz, then we clamp it to be equal to 
Csmz- This condition represents a limit on the mass ef- 
flux pvz from the disk. For sub-slow magnetosonic outflow 
from the disk Vz < Csmzi we have four boundary condi- 
tions, whereas when Vz = Csmz we have five boundary 
conditions. 

For the outer boundaries, we first consider the case 
where these surfaces are perfect conductors. Secondly, 
we consider the case of "free" outer boundaries where 
dFj/dn = on all scalar variables except for the toroidal 
magnetic field. For this field component we take [Bp • 
V(ri?0)] = on the outer boundaries, which was shown 
by Ustyugova et al. (1999) to avoid artificial coUimation 
which can come from using the "free" boundary condition 
on rB^. The free outer boundary conditions allow matter 
and Poynting fiux to freely fiow out through these surfaces. 

For the cases we discuss, the strength of the poloidal 
magnetic field at the inner radius of the disk corresponds 
to \vAp/vK)i = 16.5 and {cs/vK)i = 1, where VAp = 
|Bp|/-\/47rp. The i— subscript indicates evaluation at the 
inner radius of the disk r = Ti on the disk surface. In the 
midplane of the disk {vap/vk)z=q is less than or much less 
than unity. Different radial profiles of Cg on the disk sur- 
face have been used with similar results, including Cs/vk = 
const and Cg = const; the density profiles on the disk sur- 
face have been obtained as in Ustyugova et al. (1999). 

Figure 9 shows the evolution of the coronal plasma 
for the case of fixed, conducting outer boundaries at 
{Rm, Zra). After about 6 rotation periods of the inner disk, 
the outgoing poloidal field coUimates along the axis, 
and the returning poloidal field is pressed outwards to the 
conducting walls. Most of the magnetic field is strongly 
field- dominated with flow speeds <C w^, where va is 
the Alfven velocity. The field configuration is similar to 
that found for the high-twist Grad-Shafranov solutions as 
shown in Figure 5. 

Figure 10 shows the long time evolution of the coro- 
nal plasma for the case where the outer boundaries are 
free (i.e., open) boimdaries. These simulations evolve 
to a quasi-stationary final state where most of the re- 
gion is strongly field-dominated. In the jet region along 
the z— axis, the poloidal field is coUimated by the B^ 
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Fig. 9. — Time evolution of dipole-like field threading the disk from the initial configuration t = (bottom panels) for 
the case of conducting outer boundaries at Rm = 10 and Zm = 10. Here, ti is the period of rotation at the inner radius 
of the disk Vi. The initial field is shown by the dashed lines in Figure 1. The left-hand panels show the poloidal field 
lines which are the same as \l/(r, z) =const lines. The right-had panels show the poloidal velocity vectors Vp. For this 
calculation a uniform 100 x 100 grid was used. For times longer than ~ SU the numerical calculations crash apparently 
because of the build up of fine scale structure in the simulation region. 



field. We find that the profiles of B^{r) and Bz{r) from 
the simulations (Ustyugova et al. 2000) agree to a good 
approximation with equation (6) of our analytic model. 
The main respect in which the simulations differ from the 
Grad-Shafranov solutions is that region close to the disk 
where the magnetic field returns to the disk is not field- 
dominated. Instead this region has a hydromagnetic out- 
flow from disk. This is predicted theoretically (Blandford 
& Payne 1982) and observed in simulations (Ustyugova et 
al. 1999) to be the case when the angle 9 between the disk 
normal and the poloidal field lines is larger than 30°. 

14. CONCLUSIONS 

An ordered magnetic field threading an accretion disk 
can give powerful outfiows or jets of matter, energy, and 
angular momentum. Most of the studies have been in 
the hydromagnetic regime where there is appreciable mass 
outfiow and find asymptotic flow speeds of the order of the 
maximum Keplerian velocity of the disk, vki- These flows 
are clearly relevant to the jets from protostellar systems 
which have flow speeds of the order of vxi- In contrast, 
observed VLBI jets in quasars and active galaxies point to 
bulk Lorentz factors F ~ 10 - much larger than the disk 
Lorentz factor. In the jets of gamma ray burst sources, F ~ 
100. The large Lorentz factors as well as the small Faraday 



rotation measures suggest that these jets are in the Foynt- 
ing flux regime. This work presents self-consistent solu- 
tions for the axisymmetric, non-relativistic plasma equilib- 
ria described by the force-free Grad-Shafranov equation. 
We find solutions with magnetically collim,a,ted Poynting 
jets where there is a continuous outflow of energy, angular 
momentum, and toroidal magnetic flux from the disk into 
the external space. This behavior contradicts the com- 
monly accepted "theorem" of Solar plasma physics that 
the motion of the footpoints of a magnetic loop structure 
leads to a stationary magnetic field configuration with zero 
power, angular momentum, and flux outflows (Aly 1984, 
1991). 

Important issues remain to be investigated - the rela- 
tivistic expansion of the head of a Poynting jets into an 
external medium and the three-dimensional kink instabil- 
ity of the jet. Further, the magnetic extraction of energy 
from a rotating black hole may be important (Blandford 
& Znajek 1977; Livio, Ogilvie, & Pringle 1999). 

We thank S.A. Colgate and J.M. Finn for valuable dis- 
cussions. One of us (RL) thanks J. Chen for for valuable 
discussions about Coronal magnetic fields. This research 
was supported in part by NASA grants NAG5-9047 and 
NAG5-9735 and NSF grant AST-9986936. MMR received 
partial support from NSF POWRE grant AST-9973366. 
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Fig. 10. — Time evolution of dipolc-like field threading the disk from the initial configuration t = (bottom panels) to 
the final quasi-stationary for the case of open outer boundaries at r = i?^ and z = Z^^. Here, to^t is the rotation period of 
the disk at the outer radius Rm of the simulation region; for the parameters used, tout ^ 200ti. The initial field is shown 
by the dashed lines in Figure 1. The left-hand panels show the poloidal field lines which are the same as ^{r, z) =const 
lines; ^' is normalized by ^'o (the maximum value of 5'(r, z)) and the spacing between lines is 0.1. The middle panels 
show the poloidal velocity vectors Vp. The right-hand panels show the constant lines of —rBipir, z) > in units of 'I'o/j'o 
and the spacing between lines is 0.1. For this calculation a 100 x 100 inhomogeneous grid was used with Ar^ and Az^ 
growing with distance r and z geometrically as Ar^ = Anq^ and Azk = Aziq'^, with q = 1.03 and An = A^i = 0.05ro 
(Ustyugova et al. 2000). 
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